Gap Renormalization in Dirty Anisotropic Superconductors: 
Implications for the Order Parameter of the Cuprates 
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" Abstract. We contrast the effects of non-magnetic impurities on the properties of 

1^ superconductors having a d^2_y2 order parameter, and a highly anisotropic s-wave 

^ (ASW) gap with the same nodal structure. The non- vanishing, impurity induced, off- 

^ diagonal self-energy in the ASW state is shown to gap out the low energy excitations 

^ present in the clean system, leading to a qualitatively different impurity response of 

the single particle density of states compared to the d^2_y2 state. We discuss how this 
g behaviour can be employed to distinguish one state from the other by an analysis of 

high-resolution angle-resolved photoemission spectra. 

O 

^ PACS Numbers: 74.20.Mn, 74.25.Nf, 74.62.Dli 



The symmetry of the order parameter (OP) in the high-r^ cuprates is the most crucial 
and selective single criterion for the determination of their pairing mechanism. A number of 
recent experiments [1,2] have provided strong evidence for an anisotropic gap, possibly of 
<i^.2_y2 symmetry. The latter is the most likely candidate for a superconducting OP in purely 
electronic models, which emphasize the importance of the spin fluctuations induced by the 
strong Coulomb interactions in the CUO2 planes [3]. However, it has to be kept in mind that 
most of these experiments probe quantities which involve an average of the modulus of the 
OP over the Fermi surface. Consequently, the d^i^yi gap with its characteristic sign change 
is practically indistinguishable from a pairing state which has similar modulus, but no sign 
change, like a strongly anisotropic s-wave (ASW) state, as recently proposed [4]. Experiments 
with a direct sensitivity to the phase of the OP are therefore of greatest importance. 

Three phase-sensitive methods for the measurement of the gap have been proposed so 
far: (a) The magnetic flux dependence of the critical current of a dc SQUID constructed by 
forming weak Unks between two s-wave superconductors and the x and y faces of a single 
crystal djp._yi superconductor should exhibit a phase shift of ti, i. e., a vanishing critical current 
for zero flux, opposed to the same configuration using only s-wave superconductors [5, 6]. 
These experiments have been performed by two independent groups [7] using YBa2Cu30y_5 
as the high- -material, both suggesting a d^i-yi gap. However, similar interference experi- 
ments by other groups lead to results which seem incompatible with the dy.i_^ symmetry [8]. 
(b) Measurements of phonon Unewidths as a function of temperature were proposed [9] as a 
possibiUty to extract the k dependence of the OP, but no experimental verification has been 
reported so far. (c) A sizeable density of midgap surface states was argued to exist in a 4^2.^2 
state, whereas they should not be present in any s-wave scenario [10]. Again, there has been 
no unambigous test for the existence of these states. 

In this letter, we propose a third possible phase sensitive method which is capable of 
distinguishing a d^2_y2 from an ASW state: it is well known [11] that, opposed to an isotropic 
s-wave superconductor, scattering by non-magnetic impurities introduces severe pair-breaking 
in a state with line nodes (such as the d^i_yi state), and leads to the appearance of a finite single 
particle density of states (DOS) at the Fermi energy e^, and a rapid suppression of T^. 

As we shall show below, this behaviour is in sharp contrast to an ASW state, even if 
its nodal structure is identical with that of a d^i_yi state: due to the additional renormaUza- 
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tion of the off-diagonal self-energy [12] (which vanishes for symmetry reasons in the d^i-yi 
state), isotropic impurity scattering in an ASW superconductor leads to a smearing of the gap 
anisotropy, and the opening of a gap in the DOS over the whole Fermi surface, even if the gap 
had line nodes in the pure case. Hence, there is a clear qualitative difference in the impurity 
response of the ASW versus the <i^2_,.2 state, which should be detectable in a careful systematic 
experimental study. We propose angle-resolved photoemission (ARPES) as a suitable probe 
for this effect. However, care must be taken in the interpretation of the spectra, as we shall 
discuss at the end of this paper. 

We treat the impurity effects in the self-consistent ^-matrix approximation [13]: The Mat- 
subara components of the single particle propagator in the presence of non-magnetic impurities 
is expressed in Nambu space as (quantities with a hat represent matrices) 

^(k, iG)„) = ^ ~ , (1) 

with = co„ + il^Q, = + Sj, and E,y^- + Z3. Here co„ = Tiiln + 1)T is the Matsubara 
frequency, T the temperature. Ay. the OP (assumed to be real), and the quasiparticle energy 
measured from e^. The self-energy L = Y^]^ ^jOj (and all other matrix quantities) is expanded 
in terms of the unit and PauU matrices Gq-.-o^ and is given by '^(ieoj = rT(i(oJ, where T is 
the ^-matrix, F = n/(;rA^o)' the impurity concentration, and Nq the DOS at in the normal 
state. Finally the ^-matrix is self-consistently determined from 

fiiaj = K„ + K„Gii(Djfii(Dj, (2) 

where G(i(oJ = (kNq)~^ YlkSi^' i^n)^ ^^'^ the normal state K-matrix is parametrized by an s- 
wave scattering phase shift 5q, = -o^/c, c = cot Sq, which measures the scattering strength 
(c = unitary limit, c > 1 Bom limit). Assuming particle-hole symmetry, it can be shown 
[13] that G2 = G3 = 0. Then the non- vanishing ^-matrix components are given by 



c^-Gl + Gj ' c^-Gl + Gi 



A self-consistent solution to Eqs.(l-3) is obtained if Ay^ satisfies the gap equation 

0)„ k' ®n + ^k' + Sk' 



where Vy^, is the pair potential. 

In the following, we consider a two-dimensional electron system with a circular Fermi 
surface, and a separable pair potential Vy.y., = yT](0)Tj(0') (0 the polar angle in the plane). We 
consider two possibilities for the angular function r](0): (i) The dependence rj^i^) - cos 20, 
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which leads to the characteristic d^2_^2 gap, A(0) = A cos 20, and (ii) ri^(^) - \1 — 0|,O < 
(j) < 7t/2, periodically continued to the interval [7t/2, lit], i. e., a sawtooth function with period 
kI2. The latter case leads to a particular ASW state with A(0) = A77^(0). This is the simplest 
form of an ASW state with the same nodal structure as the d^2_y2 OP, and would therefore be 
practically indistinguishable from it in all experiments which do not probe the phase of the OP 
explicitly. In reality, it is very unlikely that an ASW OP would have real line nodes, since there 
is no symmetry reason for it to vanish. However, the existence of a small isotropic part in the 
gap would not alter the conclusions of this paper. The form of 77^,(0) was chosen to emphasize 
the different effects of impurities in the two cases, and for its mathematical simpUcity, which 
allows for the solution of Eqs. (1-4), at all temperatures and impurity concentiations. 

The fundamental difference between the two states appears in the presence of non-magnetic 
impurities, signaled by the -component of the integrated single particle propagator 



The symmetry of the d^2_y2 gap makes the sum over k vanish, and therefore Gj = Ej = 0, /. e., 
A|5^ = A,5^, whereas in the ASW state, Gj and Zj are finite. This leads to an isotropic (due to the 
assumption of pure 5-wave scattering) renormalization of the off-diagonal self-energy by the 
impurities, IS^ = A|j+Ej . Below, we shall discuss how this difference in the gap renormaUzation 
affects several observable quantities, and how this could be used to distinguish the two gap 
symmetries. 

We have self-consistently solved Eqs. (1-4) at all temperatures for both OP symmetries, 
and compared two values for the scattering strength c = 0, 1, representative of stiong and 
intermediate scattering. First consider the suppression: In both cases, we obtain the standard 
Abrikosov-Gorkov resuh [14], ln(r,o/r,) = a[^(l/2+a/(2;rr,))-^(l/2)], where a = r/(cVl) 
is the pair-breaking parameter. In the d^2_y2 case, a = 1, and we obtain the critical concentration 
= Tt^ic^ + l)NQT^/{2e^) at which vanishes, whereas in the ASW case, a = 1/4, and 
= oo. In Fig. 1(a), we show the dependence of on F (all energies are scaled in terms of 




(5) 
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Figure 1: (a) The TV-suppression as a function of impurity concentration n, for the d^-^-y^ and the ASW 
state at two scattering strengths c = 0,1. The temperature dependence of the dx2-y20P for different n„ at 
c = (b), and c = 1 (c). (d) The temperature dependence of the ASW OP for different n„ at c = 0. 



the maximum gap value Ago = A(0 -0,T - 0, rij - 0)). The value of A^q ^^ch Cu02 plane 
can be estimated to be of the order of 1 - 2 eV~^ per Cu, which leads to the relation F = vn^, 
where v ~2 meV, and n,- is in percent. The initial decrease in is given by T^-T^ = aanlA, 
i.e. the suppression in a d^_yi superconductor is four times faster than for our particular choice 
of ASW state. But note that in the ASW scenario, the coefficient a is non-generic. 

Figs, l(b-d) show the temperature dependence of the OP for different impurity concen- 
trations. In the d^2_y2 case, Figs. l(b,c), the suppression of the zero temperature gap is very 
similar to the suppression, whereas in the ASW case. Fig. 1(d), it is suppressed slower than 
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Figure 2: A comparison of the DOS at T = for the dx-^-y^ OP and the ASW OP for two scattering 
strengths (c = 0, 1) at different values of the impurity concentration. 



Next we consider the DOS, which is given by Nico)/NQ - -^tnGo(JG)„)|,-^^^^+,g. Fig. 2 
shows the impurity dependence of the DOS for the d^2_y2 and the ASW state for the two 
scattering strengths c = 0, 1 at T = 0, Note that the frequency in these plots is scaled in terms 
of Aqo, i- e., the gap suppression is taken into account. In the pure case, Nico) is linear at small 
frequencies for both states. The logarithmic singularity at <o = Aqq in the d^2_y2 case is replaced 
by a cusp in the ASW state, due to our choice of 17^(0). As soon as impurities are doped into the 
system, the difference between the two states becomes apparent: Whereas, in the <i^2_^2 state, 
the impurities lead to excitations at zero energy over the whole Fermi surface (as verified by a 
calculation of the angle resolved DOS, Ni^, (o)/Nq = 2fm{/a3„[47r^(a3,^+A^)]"^^^)}|,-^^^+,g), their 
effect in the ASW state is to induce a finite energy gap, which initially grows substantially as 
a function of the impurity concentration n,. For very large n,, the gap anisotropy in the ASW 
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state becomes completely smeared out by the isotropic scattering, and the DOS approaches the 
familiar shape of an isotropic s-wave gap. This should also have the consequence of activated 
behaviour in thermodynamic and transport properties of dirty samples at low temperatures, in 
contrast to the power laws predicted for clean [15] or dirty [16] d^i_y2 states. Note also that 
the resonance behaviour found in the unitary scattering limit (c = 0) of the <i^2_^,2 state at low- 
frequency is cut off in the ASW state. This is due to the fact that in the latter case, the finite gap 
renormaUzation prevents the denominator of Tq in Eq. 3 to become singular. Consequently 
there is no quaUtative difference between the results obtained in the Bom (c > 1) and unitary 
scattering limit, and all the results can essentially be parametrized by the effective normal state 
inverse scattering rate r/(c +1). 

The indication of an anisotropic gap found by ARPES [2] motivated us to calculate the 
impurity dependence of the spectral function A(k, (O), which is defined as 

A(k, 0)) = —3m {gii(k, ia)r)\ico„=co+i8] - (6) 

where is the 11 component of the single particle matrix propagator. In the pure case of the 
BCS models we consider, the spectral function consists of the familiar two 5 functions at the 
quasiparticle energies (O = ±Ey^ [17] 

A(k, (0) = uldico - £k) + vidiOD + E^.), (7) 

where Ey^ - ^J^i. + A^, = (1 + (^^/£'^)/2, = (1 - ^-JE^Il. Hence the quasiparticle peaks 
shift by an energy oiE^-\Eiy\ when going from the normal to the superconducting state. 

Before we proceed with the results obtained in the presence of impurities, let us briefly 
discuss the relation of the spectral function to the intensity /(k, ©) measured in an ARPES 
experiment. For a perfectly flat surface, the photoexcited electron preserves its momentum 
parallel to the surface, when leaving the sample. In the case of the cuprate superconductors, 
one therefore cleaves the crystal with a surface parallel to the a - b plane, so that momen- 
tum information is obtained about the relevant CuOj planes. Then the intensity becomes a 
function of the normal quasiparticle energy (§, the angle 0, and (O. The energy difference of 
the outgoing electron to the incoming photon gives information about the negative frequency 
part of A(0, ^, CO). The substantial background intensity which is usually observed is caused 
by multiple inelastic scattering of some photoexcited electrons before they leave the sample. 
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Figure 3: The averaged (bottom curves) and true (top curves) spectral functions for the dx2-y2 and the 
ASW OP at 1^ = 0, = for varying impurity concentration in the unitary limit. 



For a comparison of theoretical results to the measured spectra, it is important to take into 
account the finite experimental resolution, both in energy and momentum. We have simulated 
the effect of finite angular resolution by a Gaussian width for the normal quasiparticle en- 
ergy, and similarly the finite energy resolution by a corresponding width A^, for the frequency. 
We then argue that the measured intensity should be proportional to A(0, ^,(o<0) integrated 
with these probabiUty distributions. 

In Fig. 3, we plot the spectra for A^kp, = 0, ct) < 0), i. e., along the direction of the 
gap maximum, calculated as a function of impurity concentration in the unitary limit for both 
states. We always plot the spectra for -0, i. e., at the Fermi surface, and T -0. The top and 
bottom curves are obtained before and after the averaging process respectively. The resolution 
widths were chosen as A^/Aqq = 0.2, and A^/Aqq = 0.5. Assuming a value of Ago ~ 30meV, 
these values correspond to the best available resolutions on current spectrometers. 
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Figure 4: The averaged (bottom curves) and true (top curves) spectral functions for the d^^-i^ and the 
ASW OP at 1^ = 0, (/» = 7r/4 for varying impurity concentration in the unitary limit. 



For increasing impurity concentration the effect on the spectra is qualitatively similar for 
both states: the quasiparticle peaks shift to smaller energy (due to the reduction in Aq), and they 
become substantially broadened. Around zero frequency, the spectra are different: Similar to 
the total DOS, the spectral weight is non-zero in the d^i_^i case, whereas there is a true gap 
in the ASW state around © = 0. However the broadened spectra indicate that this difference 
would be difficult to detect, especially since it occurs in the vicinity of e^, where the ARPES 
experiments suffer from a sharp drop in intensity. 

Fig. 4 shows the spectra along the node directions calculated with the same parameters 
as in Fig. 3. Along this direction, the two states show a clear qualitative difference: for the 
<i^.2_,.2 state, the only effect of an enhanced impurity concentration is the broadening of the peak 
which remains right at e^. In the ASW case, however, the opening up of a gap in the presence 
of impurities makes the peak shift away from e^. When looking at the averaged spectra, we 
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notice that the shift of the peak is still clearly visible, even though its sharp cutoff is smoothed 
out by the finite resolution. This shows that a clear distinction of the ASW from the d^i^^i state 
should be detectable in a systematic analysis of high-resolution ARPES spectra taken at the 
possible line nodes as a function of non-magnetic impurity concentration. Note that a similar 
effect would happen if the OP in the ASW state had a small isotropic part. 

Since the experiments are done at finite temperature, it is important to divide the spectra 
by the Fermi occupation factor (after a careful background subtraction) for a comparison of 
the data to our theoretical curves. Figs, l(b-d) show that the reduction of A by increasing the 
temperature is negligible if TIT^. < 0.2 (as in actual measurements), so that our T = results 
should be applicable when thermal occupation of quasiparticle states with £"1^ > is taken into 
account. 

A final comment concerns the data published by Shen et al. [2]. They showed that the 
anisotropy of their extracted gap values decreases as a function of time after the cleavage of 
the sample. If this aging would be due to the creation of non-magnetic impurities, e. g., by 
the diffusion of oxygen, then the effect could be understood in terms of an ASW state as we 
showed above. It is also possible that it is due to a degradation of the surface quality, as Shen 
et al. argue, in which case the resulting poorer k resolution could be responsible for the effect. 
But in this case, one might expect (i) an enhanced width of the superconducting quasiparticle 
peak which was not observed, and (ii) that the poorer k-resolution should also wash out the 
k-dispersion in the normal state data. A sytematic study of bulk irradiated samples would 
certainly be more conclusive, and is therefore highly desirable. However it has to be made 
sure that the radiation damage does not lead to localized magnetic moments, or a shift in the 
chemical potential, since this would mask the effects we were discussing. 

In conclusion, we have compared the effect of non-magnetic impurities on the properties of 
superconductors with a d^2_y2 versus a highly anisotropic s-wave OP, both with the same nodal 
structure. We have shown that even though the two states are indistinguishable in experiments 
which probe Fermi surface averages of the OP, their behaviour is qualitatively different in the 
presence of non-magnetic impurities. The difference comes from the non-vanishing, impurity 
induced, off-diagonal self-energy in the ASW state, which leads to the opening up of a true 
gap in the single particle excitations of dirty samples. In particular, we discussed how this 
effect could be observed in ARPES experiments. 
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